## This code creates Figure 4.

## This code requires Google Mpas API key to run. Please run register_google([Insert Google API Key]).

mm_basemap <- get_googlemap("Phyu, Myanmar", zoom = 7, maptype = "terrain")

# Figure 4

ggmap(mm_basemap) +
  geom_polygon(aes(x=long, y=lat, group=group), fill='white', size=1,fill='gray', data=fortify(bago),alpha=0.6) +
  geom_point(data = as.data.frame(colonial_towns_shp[colonial_towns_shp$replaced==0,]),aes(x = longitude, y = latitude),colour = "black",shape=4,size=2) +
  geom_point(data = as.data.frame(colonial_towns_shp[colonial_towns_shp$replaced==1,]),aes(x = longitude, y = latitude),colour = "black",shape=16,size=2)
